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ABSTRACT 



Context. Stars in young nearby associations are the only targets allowing giant planet searches at all separations in the near future, 
by coupling indirect techniques such as radial velocity and deep imaging. These stars are first priorities targets for the forthcoming 
planet imagers on 8 to 10 meter class telescopes. Young stars rotate more rapidly and are more active than their older counterparts. 

^vi Both effects can limit the capability to detect planets using RV. 

Aims. We wish to explore the planet detection capabilities of a representative sample of stars in close and young associations with 
radial velocity data and explore the complementarity between this technique and direct imaging. 

i i Methods. We observed 26 such targets with spectral types from A to K and ages from 8 to 300 Myr with HARPS. We compute the 

detection limits with two methods, in particular, a method we recently developped, that takes into account the frequency distribution 
of the RV variations. We also attempt to improve the detection limits in a few cases by correcting for the stellar activity. 
Results. Our A-type stars RV show high frequency variations due to pulsations, while our F-K stars clearly show activity with more 
or less complex patterns. For F-K stars, the RV jitter and v sini rapidly decrease with star age. The data allow us to search for planets 
ft with periods typically ranging from 1 day to 100 days, and up to more than 500 days in a few cases. Within the present detection 

~. limits, no planet was found in our sample. For the bulk of our F-K stars, the detection limits fall down to sub-Jupiter masses. We show 

that these limits can be significantly improved by correcting even partially for stellar activity, down to a few Neptune masses for the 
least active stars. The detection limits on A-type stars can be significantly improved, down to a few Jupiter mass planets, provided an 
appropriate observing strategy. We finally show the tremendous potential of coupling RV and AO deep imaging results. 
Conclusions. The RV technique allows the detection of planets lighter than Jupiter, and down to a few Neptune masses around young 
stars aged typically 30 Myr or more. Detection limits increase at younger ages, but (sub-)Jupiter mass planets are still detectable. In 
the forthcoming years, the use of complementary techniques will allow a full exploration of the Jupiter mass planets content of many 
^- of these stars. 
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• , 1. Introduction diate period planets. This suggests that dynamics, due for exam- 

. . pie to planet-planet interactions or early disk-planet interactions, 

Thanks to hundreds of discoveries (http://exoplanet.eu/) for plays an important role in the building of planetary systems. 
more than 15 years, mainly coming from radial velocity (RV) 

and transit surveys, our knowledge of exoplanets has dramati- Giant planets play an important role as they carry most of 

H cally improved. We know that exoplanets are frequent around the planetary system mass. They therefore strongly impact the 

solar-type, main sequence (MS) stars: more than 50% have plan- dynamics and fate of lighter planets and the final architecture of 

ets (all masses, Mayor et al., 201 1) and about 15% have planets the planetary systems. Giant planets also impact the detectabil- 

with masses larger than 50 M Earth . An unexpected diversity of Hy of lighter planets as far as indirect methods are concerned, 

planet properties such as separations, eccentricities and orbital Even though giant planets represent the majority of the planets 

motions (for example retrograde orbits discovered thanks to the detected so far, we are far from having a clear picture of their oc- 

Rossiter-McLaughlin effect) was revealed for short or interme- curence, variety and properties. RV and transit explorations are 

indeed still limited to planets orbiting within 5 AU from their 
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* Based on observations made with the ES03.6m/Harps spectrograph P arent stars > whlle current deep (adaptive-optics, herefater AO) 

at La Silla.083.C-0794(ABCD); 084.C-1039(A); 084.C-1024(A) imaging, meant to detect more distant planets, are not sensitive 

** Tables of radial velocities are only available in electronic form at enough for MS stars, except around a few young early-type stars 

the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via which have already reached the MS. Stars in young nearby as- 

http://cdsweb.u-strabg.fr/cgi-bin/qcat?J/A+A/ sociations (see for example Torres et al., 2008, Zuckerman & 
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Song, 2004) are therefore the best targets for a complete giant 
planet exploration, combining the data obtained with forthcom- 
ing planet imagers such as SPHERE at the VLT (Beuzit et al., 
2008) or GPI at GEMINI (Macintosh et al., 2008) and high pre- 
cision spectrographs. In addition, these data could be combined 
in the future with astro me trie measurements when available. 

Compared to mature, MS stars, limited effort has been de- 
voted so far to the search for planets around young stars with RV 
techniques. To our knowledge, the only significant surveys were 
performed by Paulson et al. (2004) on 94 stars members of the 
Hyades (aged about 600 Myr) and on 6 1 stars members of var- 
ious moving groups (MGs) or stellar associations aged between 
12 and 300 Myr (Paulson & Yelda, 2006). In both cases, high 
spectral resolution data, in the range 40000-70000 were used. 
The internal errors were 3-5 and ^ 10 m/s. The Hyades stars 
were found to be weakly active, with stellar jitters in the range 
10-20 m/s (average value = 16 m/s). No planet was detected with 
masses down to 1-2 MJup, and down to the sub-Jupiter mass 
regime in a few cases and with periods typically of 1-2 years 
or less. In the second case, no planet was detected down to 1-2 
MJup, for periods up to 6 days for the oldest stars while the de- 
tection limits fell into the brown dwarf regime for the youngest 
targets (members of the /? Pic group). 

Stellar activity can mimic planet signatures (Desort et al., 
2007), and great care must be taken when analysing RV data. 
In the past, two detections were announced around TW Hya 
(Setiawan et al., 2007), and BD +20 1790 (Hernan-Obispo 
et al., 2010), but were strongly debated (Figueira et al., 2010). 
Bailey et al. (2012) used NIRSPEC at Keck to observe a set 
of 20 young stars with high precision spectroscopy in the near- 
IR, where line shifts due to spot-induced activity are expected 
to be smaller than at optical wavelengths. Detection limits of 
8 Mj up (respectively 17 Mj up ) were obtained for periods of 3 
days (respectively 30 days) periods. Crockett et al. (2012) used 
CSHELL for similar purposes on 9 T Tauri stars and showed that 
for these stars, the jitter in the near IR was reduced by a factor 2- 
3 compared to the jitter measured in the optical range. No planet 
was firmly detected. Finally, we detected unambiguously a 2.8 
Mj up giant planet orbiting with a 320 days period, a 150 Myr, 
weakly active, star (Borgniet et al., 2013). 

We investigate here the potential of high precision optical 
spectroscopy to search for giant planets around 26 stars which 
are members of close-by young associations, as such stars are 
potential targets of the SPHERE NIRSUR survey aimed at imag- 
ing giant planets. Section 2 describes our targets and observing 
log. Section 3 presents the data analysis. The results are shown 
and discussed in Sect. 4. In Section 5, we show how such re- 
sults can be coupled to deep imaging detection limits to perform 
a complete exploration of giant planets in the environements of 
these young and close stars. 



2. Targets and observations 

2.1. Target list 

The bulk of our sample is made of 22 bright (V < 9), nearby 
(d < 80 pc), young (< 200 Myr) stars identified in young 
MGs and from systematic spectroscopic surveys. Four other 
stars with ages greater than 200 Myr were observed: HD90905, 
HD109536, HD105690 and HD188828. The 26 stars are listed 
in Table 1. Our sample includes 8 A-type stars, 4 F-type stars, 
10 G-type stars and 4 K-type stars. Several of them are known 
to have debris disks (Rodriguez & Zuckerman, 2012). In sev- 



eral cases, peculiar disks properties (shapes, offsets) suggest that 
inner planets could be present. 

We provide now some details about a few peculiar targets. 
Otherwise indicated, relative disk luminosities Ldi s k/L star are 
from Moor et al. (2006) or Rhee et al. (2007), and binary in- 
dications from Rodriguez & Zuckerman (2012): 

- HD39060 (J3 Pictoris ) is surrounded by a seen edge-on disk 
(Ldisk/L star = 2.10~ 3 ), as well as a planet detected in direct 
imaging (Lagrange et al., 2009b, 2010). Detection limits in 
deep imaging fall in the giant planet regime. A study of RV 
data allowed Lagrange et al. (2012) to set an upper limit to 
the dynamical mass of /3 Pictoris b. 

- HD61005 (The Moth) is surrounded by a thin and narrow 
disk which is resolved in near-IR (Hines et al., 2007, Buenzli 
et al., 2010, and references therein). The disk is seen as a ring 
inclined by 84 degrees with respect to pole-on and is offset 
from the star center. In addition, streamers are observed at 
the edge of the ring. Detection limits in deep imaging fall in 
the giant planet regime. 

- HD172555 belongs to the (3 Pictoris MG (Zuckerman et al., 
2001). It has a spectral type similar to that of /3 Pictoris but an 
IR excess smaller than that of /? Pictoris (Ldisk/L sta r= 8.10~ 4 ). 
It is also part of a wide binary system. 

- HD174429 (PZ Tel) has an IR excess (L disk /L star = 9.10~ 5 ). A 
brown dwarf companion with a 20-40 Mj up mass is present 
with a projected separation of 0.3" representing 15 AU 
(Biller et al., 2010). Given the companion AJ (5.6 mag), and 
its expected V-J (5.5-6.5 mag) given its mass and age, and 
the V-K ( 1 .2 mag) of the star itself, the contrast at optical 
wavelengths between PZ Tel and its companion is about 9-10 
mag. The companion contribution to the visible flux within 
the Harps fiber (1") is therefore negligeable. 

- HD 18 1327 is a member of the f3 Pictoris MG, with a large IR 
excess (Ldi s k/L star = 2.10~ 3 Lebreton et al., 2012). A debris 
disk has been resolved in near IR by Schneider et al. (2006) 
as a ring located at 90 AU and inclined by about 30 degrees 
from face -on. 

- HD2 16956 (Fomalhaut) is also surrounded by a debris disk 
(L,ji S k/L star = 8.10~ 5 ) seen as a narrow ring at 1 15 AU (Kalas 
et al., 2005), inclined by 24 degrees from edge-on and offset 
from the star. It may also be surrounded by a planet orbiting 
just inside the ring (Kalas et al., 2008), and shaping its sharp 
inner edge (Chiang et al., 2009), but this planet is still de- 
bated (Janson et al., 2012). Recently Lebreton et al. (2013) 
showed that warm dust was also present within AUs from the 
star. 

- HD218396 (HR8799) is surrounded by a system of 4 planets 
detected in imaging with projected separations between 15 
and 70 AU (Marois et al., 2008, 2010). A thin debris disk 
(Ldi s k/L st ar= 2.10~ 4 ) has also been resolved at sub-mm wave- 
lengths (Patience et al., 201 1) and has an inclination of 20-50 
degrees from face-on (Kalas et al., 2010). 

2.2. Observations 

We have obtained more than 2000 high resolution (R ^ 1 15000) 
and high signal to noise (S/N) spectra of our targets with the 
fiber-fed spectrograph Harps (Mayor et al., 2003) at the La Silla 
3.6m telescope. In the case of fi Pictoris, we used in addition to 
our data, a few data available in the archive and obtained with 
the same set-up. The spectra cover a wavelength range between 
3800 and 6900 A. During our runs, we usually recorded two 
consecutive spectra for each telescope pointing. In the follow- 



A.-M. Lagrange et al.: Planets around stars in young nearby associations 



ing, one pointing will be refered to as one "epoch". As our main 
goal was to look for short period planets (100 days or less), we 
recorded several spectra of each target during each run whenever 
possible. In addition, we also recorded whenever possible spec- 
tra at two or three different epochs during one night, to identify 
possible high-frequency RV variations. For the brightest stars, 
we did not consider spectra with saturation that could affect the 
RV measurements. We report in Table 2 information about each 
star observation. The time-span varies between a few hours for 3 
of our targets and 1912 days. The stars with a long time-span are 
early-type stars which happened to be part of a previous survey 
on A-F stars described in Lagrange et al. (2009a). The exposure 
times were computed to get a S/N larger than 200 at 550 nm 
in most cases. This corresponds to exposure times between less 
than 1 minute for the brightest targets and 15 minutes for the 
faintest ones. 

Finally, for some of the earliest types, generally pulsating 
stars, we performed continuous observations over typically 1-2 
hours to check that short-term pulsations are indeed present and 
to estimate their amplitude. 



3. Data analysis 

3.1. RV measurements 

The RV were computed using SAFIR, a tool dedicated to mea- 
sure accurate RVs on fast rotators, described in Galland et al. 
(2005), and based on the Fourier interspectrum method devel- 
oped by Chelli (2000). SAFIR has been extensively used, in par- 
ticular in RV surveys of rapidly rotating stars (Lagrange et al., 
2009a). SAFIR input data are the extracted 2D spectra provided 
by the Data Reduction Software (DRS) pipeline, and the re- 
lated wavelength calibration coefficients. SAFIR takes care of 
the blaze correction as well as bad-pixel correction (Galland 
et al., 2005). The errors associated to the individual RV mea- 
surements range between typically 1 .5 m/s for the slowly rotat- 
ing stars to typically 30-40 m/s for the fast rotating stars. The 
values obtained for slowly rotating stars are similar to those ob- 
tained by the HARPS DRS; these values show that the instru- 
ment internal errors are very low. Lomb-Scargle periodograms 
of the RV variations were also computed for each target. 

3.2. Line profile variations and activity indicators 

For stars with v sin; below typically 100 km/s and spectral types 
later than A3V, it is possible to compute the bisector velocity 
spans (BVS) using the cross-correlation function, as shown in 
Lagrange et al. (2009a). Typical examples are provided in Fig. 1 
for an active star (upper panels) and a pulsating one (lower pan- 
els). The relation between BVS and RV is useful to further 
characterize the temporal RV variations whenever present. A 
correlation between the RV and BVS variations indicates that 
the RV are due to stellar activity (Desort et al., 2007) while 
pulsation-dominated variations induce instead uncorrected (RV, 
BVS) variations. We also computed the v sin; using the cross- 
correlation function for the low to moderately rotating stars and 
the obtained values are given in Table 2. 

For active stars, we also computed the S-index that measures 
the amount of chromospheric activity, and derived the Ca index 
R'hk according to Duncan et al. (1991) 1 . In the case of the Sun, 



there is a correlation between the long-term Ca emission vari- 
ability and the RV variability (see Meunier & Lagrange, 2013, 
and references therein). For solar-type MS stars, a correlation 
between the long-term photometric variability and the Ca emis- 
sion is also observed (see in particular Lockwood et al., 2007): 
the brightness increases when the Ca emission increases. This 
is consistent with the fact that the long-term photometric varia- 
tions are, as for the Sun, dominated by plages. Correlations be- 
tween RV and Ca emission have been searched for (Santos et al., 
2010, Isaacson & Fischer, 2010, Boisse et al., 2009, Dumusque 
et al., 2011), but only rarely reported (Dumusque et al., 2011, 
and references therein). However, instrumental artefacts such as 
diffused light could alter the Ca emission measurements and/or 
the RV measurements (the long term stability must be ensured 
to a few m/s) and may therefore prevent finding such a correla- 
tion. For young solar-type stars, the situation is even less clear, 
as much less data are available. The long-term photometric vari- 
ations are anti-correlated with chromospheric activity variations, 
i.e. the brightness decreases when the Ca emission increases, in 
contrast to older Sun-like stars (Lockwood et al., 2007). This is 
interpreted as the photometric variations being spot-dominated. 
For stars with a high enough v sin; (i.e. larger than 15-20 
km/s) or with very strong emissions, the 1 A wavelength range 
used to measure the Ca emission index is too narrow compared 
to the emission itself. Even though sophisticated modeling has 
been proposed to estimate the Ca index in such cases (Schroder 
et al., 2009), we did not make any attempt to do so, as the main 
purpose of the present paper is not the study of stellar activity. 
In such cases no value is given in Table 2. 

3.3. Detection limits measurements 

As extensively discussed in a previous dedicated paper (Meunier 
et al., 2012, hereafter MLB 12), there are different ways to es- 
timate detection limits associated to RV time series. We will 
use two methods in the present paper, namely the "rms-based 
method" and the LPA method. Noticeably, we will not use the 
bootstrap method which, we showed, provides detection limits 
not as good as the LPA method and in the case of pulsating stars, 
may provide unrealistic results. We refer to MBL12 for a de- 
tailed description of these methods and on their respective mer- 
its. We briefly summarize the principles of the "rms-based" and 
LPA methods. Examples of obtained detection limits are pro- 
vided in Fig. 2. 

The fastest method, the "rms-based method" is based on the 
measured standard deviation (hereatfer rms) of the RV serie and 
was first described in Lagrange et al. (2009a). In the rms-based 
method, we compute the RV produced by a planet with a given 
mass and period for the actual temporal sampling and then the 
rms of these simulated RV. This is done for 1000 orbital phases. 
The distribution of this rms is Gaussian. If the rms of the ob- 
served RV is lower than the average rms of the 1000 simulated 
RV, the planet is detectable. The level of confidence (detection 
probability) is obtained by comparing the standard deviation of 
the simulated distribution and the difference between the ob- 
served rms and the average value of the simulated distribution. 



1 Note that we did not apply any correction factor to the obtained 
K'hk index, as it is sometimes done to calibrate the values with those 
of Mount Wilson, because 1) we are mainly interested here in relative 



variations, and 2) such corrections are possible only when a lot of values 
(for several targets, with stable levels of activity), taken both at Mount 
Wilson and in the other observing site are available. Here, we are pos- 
sibly dealing with a strong variability, so a correction could be applied 
only if we had simulatenous measurements, for the same targets. Our 
log(R'# K ) are similar within 0.3 dex with those already published ones 
however. 
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Fig. 1. Upper panels: RV versus time, bisectors, BVS versus time and BVS versus RV for the active G-type star HD 1834 14. In the last panel, we 
indicate the slope of the (RV, BVS) correlation. Lower panels: same for the pulsating A-type star HD216956. 



We showed in Lagrange et al. (2009a) that provided a very good 
temporal sampling, the rms-based detection limit in mass for a 
given period is similar to the mass of a planet that would produce 
an RV amplitude equal to 3xrms(RV). In the following, the later 
limit will be referred to as the rms-based achievable limit. 

The LPA (local power analysis) method provides more ac- 
curate and lower detection limits when enough data are avail- 
able (see below). Briefly, for a given period P, we compute the 
maximum of the periodogram of the observed RV in the range 
0.75P-1.25P 2 : this provides a threshold for that period. Then, for 
a given planet (mass, period), we compare the maximum ampli- 
tude of the power spectrum of the RV induced by this planet on 
the same temporal sampling with that threshold, for 100 trials 
of the phase of the planet. If for all these trials, all "planet am- 
plitudes" are above the threshold then the planet mass is above 
the detection limit. We iterate on the planet mass until we reach 
a mass for which this condition is no longer met, indicating a 
mass below the detection limit. The smallest step we use is in 
general 0. 1 Mj up , but we chose a smaller step in a few cases to 
avoid a discretization of the detection limits. With our present 
set of data, we can use the LPA method for ten stars, for which 
more than 40 data points were available. 

As expected, the LPA detection limits are much lower than 
the rms-based limits. To quantify the improvement, we compute 
the ratio between the rms-based and the LPA detection limits 
for each of the 200 periods considered, and then the median 
of the 200 ratio. The median ratio is typically a factor of 2 to 
3, but may sometimes reach higher values (see Section 4). The 
improvement is due to the fact that the LPA method takes into 
account the variability of the power amplitude frequency, either 
due to an actual signal or to the temporal sampling, while the 
rms-based does not. The advantage of the rms-based is that the 
computation is fast, and it is less sensitive to the number of data 
available. It provides an order of scale but pessimistic estimation 
of the achievable detection limits (Lagrange et al., 2009a). 

3.4. Search for circumstellar Ca in rapid rotators 

Finally, as a by-product of these observations, we also used the 
spectra to search for narrow Call absorption lines at the bot- 
tom of the rotationnally broadened stellar lines. Such absorp- 



tions could indicate the presence of circumstellar gas, as present 
in the case of HD39060 (e.g. Lagrange et al., 2000). 



4. Results 

4.1. RV time series and variability properties 

We show in Appendix A the obtained RV times series for our 
26 targets. The RV rms and amplitudes, as well as the BVS rms 
and amplitudes for stars with enough data points, are reported in 
Table 2. 

The average jitter of our eight A type stars is 237 m/s, while 
that of the three F-type stars is 35 m/s and that of the thirteen G- 
K stars is 13 1 m/s. The median is 67 m/s in the latter case. For the 
stars for which the BVS could be computed, the average BVS is 
61 m/s for the F-type stars and 74 m/s for the G-K type stars. We 
recall though that within a given spectral type range, there is a 
large dispersion of jitters due to the star ages (see below). Among 
the 14 F-K stars with at least two epochs 3 , only one, HD224228, 
shows no detetectable level of RV variations, with a RV jitter 
below 5 m/s. Twelve stars have jitters larger than 5 m/s. Among 
them, only three have jitters larger than 150 m/s (they will be 
refered to as "high amplitude variable stars"), five have jitters 
between 50 and 105 m/s (they will be refered to as "medium 
amplitude variable stars") and four have jitters between 5 and 50 
m/s (they will be refered to as "low amplitude variable stars"). 

The v sin/ deduced from our cross-correlation functions are 
also provided in Table 2, whenever possible. For the stars 
with already published values, our values agree within typi- 
cally 10 % with previously reported ones, except for the low 
v sin/ HD224228 for which we find 5.7 km/s instead of 8.7 km/s. 
The average v sin/ is 20. 1 km/s for F-type stars and 20.5 km/s for 
G-K stars (with a median of 13.3 km/s). 

We also provide, whenever possible, the R'hk values. They 
agree within 0.1-0.2 dex with reported ones, although the lat- 
ter often show discrepancies from one author to the other. The 
differences can be due either to systematics or to intrinsic stel- 
lar variability, which would not be surprising as all these stars 
(except HD224228) are active. 

Finally, in Fig. 3, we show the RV rms and v sin/ as a function 
of the star ages, as well as the v sin/ as a function of RV jitters. 



2 We checked that considering different period ranges, such as 0.7P- 
1.3P, 0.8P-1.2P, 0.9P-1.1P does not significantly change the detection 
limits curves. 



We do not consider HD 177 171 in this section dedicated to the ori- 
gin of the RV jitter, but we will consider it in the next section to get 
estimates of detection limits on this type of stars as well. 
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Fig. 2. Upper panel: RV versus time, periodogram of RV and detection limits [from left torighi) for the G2V star HD183414, coverin: 
of 189 days. The detection limits have been computed using the rms-based method (black line), LPA (orange stars). Lower panels: 
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1000 



; a time span 
same for the 



500.0 ■ 


■ 














? 400.0 ■ 

g 300.0 ■ 

E 

200.0 ■ 


♦ 




■ 










♦ 

I 

♦ 














100.0 ■ 
0.0 ■ 


■ 
♦ 


♦ 

♦ 
-* — , — 




♦ 


♦ 


■ 


A 



100 



200 300 

Age (Myr) 




200 300 

Age (Myr) 



500 



being a short period spectroscopic binary, with a possible period 
of 1.7 days (we note that HD 177 171 was classified as a close 
binary by Frankowski et al., 2007, on the basis of Hipparcos as- 
trometric data). We see that globally, the v sin/ of late-type stars 
significantly decreases with age. This is in qualitative agreement 
with the results of Weise et al. (2010) analysis on a larger sam- 
ple of stars and covering a larger age range. The RV jitter also 
globally decreases with age, and increases with v sin/ . For early- 
type stars, the correlations between v sin/ and age, or RV jitter 
and age, are not as clear as for late-type stars. 



400 500 4.2. Origin of the variability 



Fig. 3. Upper panel: jitter versus age for our stars with enough data 
points (except HD177171), for G-K stars (blue diamonds), F-type stars 
(green triangles) and A-type stars (red squares). Note that an age of 300 
Myr was assumed for Fomalhaut. Bottom panel: same for v sin; versus 
age. 



We considered only stars with enough data available. We did not 
consider HD 1 77 1 7 1 (F6 V) either, because it shows clear signs of 



We used the bisectors (see Appendix A) as well as extracted 
short-term time series of pulsating stars (see examples in Fig. 4) 
to classify the observed variations. The results of the classifica- 
tion (activity, pulsations, spectroscopic binarity) are reported in 
Table 2. 



4.2.1. Early-type stars 

All our A-type stars show high frequency and high amplitude (up 
to 5 km/s) RV variations which are characteristic of pulsations 
(see Table 2). Similar conclusions were reached in Lagrange 
et al. (2009a) for MS stars with similar spectral types. 



4.2.2. Late-type stars 

The RV and BVS of late-type stars are generally well corre- 
lated, as shown for example for HD183414 in Fig. 1, indicating 
that the variability is due to stellar activity. We do not find any 
correlation between R'hk and RV for any of our stars. This is 
compatible with the long-term activity of late-type stars beeing 
dominated by spots instead of plages, as proposed by Lockwood 
et al. (2007). Our observations show a similar effect on shorter 
timescales. 
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Fig. 4. Upper panel: High frequency variations of RV for HD39060 (ex- 
posure time of 90 seconds). Middle panel: same for HD216956 (expo- 
sure time of 8 minutes). Lower panel: same for HD218396 (exposure 
time of 40 seconds). 



To help interpreting our results, we ran several cold spots 
simulations using S AFIR (see a description of the simulations in 
Desort et al., 2007) to characterize the RV and BVS variations 
using simulated spectra over the 377-691 nm region. All these 
simulations are done assuming a solar-type star, but we checked 
that the results for F5 and G2 stars are quite similar. The spot 
temperature contrast was 1200 K (Berdyugina, 2005). Also, in 
these examples, a very high S/N (> 300) was considered, to fo- 
cus on the effects of star and spot parameters. Results are pro- 
vided in Fig. B. 1 in the case of a star seen edge-on with different 
vsin/, and in Fig. B.2 for different star inclinations, vsin/ and 
spot latitudes. 

In most cases (HD37572, HD45270, HD90905, HD102458, 
HD141493, HD174429, HD181327, HD183414, and 
HD2 17343), the (RV, BVS) diagrams indicate activity with a 
simple pattern, characteristic of either a single spot or slowly 
evolving spots at the same latitude (Desort et al., 2007). A more 
complex case, HD 105690, is described below. 

We detail here a few cases of interest. 



to determine any periodicity in these variations, but they are 
compatible with a 5 day period, as deduced from photomet- 
ric data by Desidera et al. (201 1). Setiawan et al. (2008) ob- 
served this star with FEROS and reported strong variations 
with an amplitude of 150 m/s. Within the 10 m/s FEROS 
precision on RV, they found a correlation between RV and 
BVS as an indication of the presence of spots. Our results 
are compatible with these conclusions. 
HD 105690: the RV and BVS data shown in Fig. 5 are rather 
puzzling and very different from those of HD 1834 14 for ex- 
ample (Fig. 1). Most of the data are compatible with a single 
spot with a period of 4.88 days, as checked by fitting the RV 
data by a set of two signals with periods 4.88 and 2.44 days, 
which leads to a correlation between RV and BVS. However, 
the data recorded between JD 2454343 and JD2454348 lead 
to a more complex (RV,BVS) diagramme, with an additional, 
more inclined and off-centered contribution. The RV varia- 
tions between JD 2454343 and JD2454348 are nonetheless 
still compatible with a period of about 4.8 days. The com- 
parison with spot simulations shows that the observed (RV, 
BVS) characteristics can be reproduced only when assum- 
ing a small star inclination, of the order of ^ 10-20 degrees 
with respect to pole-on. If we do not consider the data taken 
between JD 2454343 and JD2454348, the observed RV are 
compatible with a spot with a projected surface of about 1% 
at =* 30 degrees from the star pole. The spot latitude is not 
well constrained as the star is close to pole-on. Obviously, 
this spot must be long-lived. This is in agreement with our 
knowledge on young stars, even though short-lived spots can 
sometimes be found as well, as shown by Garcia-Alvarez 
etal. (2011). 

We tried to see whether we could reproduce the RV and BVS 
properties for the whole period by adding a short-lived spot. 
We did not find any configuration that would produce an ac- 
ceptable match. We can therefore exclude additional spots 
at a similar latitude as the main, long-lasting spot, or a spot 
at a lower latitude. More data and studies will be needed to 
understand HD 105690 RV variations. 
HD 174429: it shows high amplitude RV and BVS variations 4 
(see Fig. 6). The observed bisector shape is closer to that 
of pulsating stars than to that of low v sin/ late-type active 
stars. Simulations show that for a given spot size, the RV and 
the BSV amplitudes significantly increase with vsin/ (see 
Fig. B.l and Fig. B.2, so we conclude that the high RV and 
BVS amplitudes can be just a consequence of the star high 
v sin/ . Note that the PZ Tel v sin/ , 80 km/s, is the highest pro- 
jected rotational velocity among our late-type stars. Garcia- 
Alvarez et al. (2011) proposed that this high vsin/ is due to 
accretion from a close TT star. We do not see spectroscopic 
signs of accretion, though. 

HD207575: its (RV, BVS) shape is quite different from that 
of our other active F stars, and closer to that of earlier type 
stars. However, for a given spot pattern, the RV and BVS 
characteristics strongly depend on vsin/, as illustrated in 
Fig. B. 1 and Fig. B.2. The (RV, BVS) shape departs from the 
roughly linear or even from the "eight-shape" form for stars 
seen with very low (5 degrees) inclinations from pole-on and 
with vsin/ in the range 30-40 km/s. Such a situation is how- 
ever not compatible with the star spectral type and observed 
vsin/ of 35 km/s. 



HD61005: only three epochs are available, spread over a 
few days only. The data reveal variability, with an ampli- 
tude larger than 100 m/s. The amount of data is not sufficient 



Note that the available data do not allow to detect the trend due to 
PZ Tel B, due to the limited time span available, and the amplitude of 
the RV variations. 
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Fig. 5. From left to right: RV variations of HD 105690, bisectors, bisector velocity span (BVS) versus time and BVS versus RV. In the (RV, BVS) 
diagram, the data corresponding to the period from JD 2454343 to JD2454348 are in red. 
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4.3. Detection limits 

4.3.1. Obtained detection limits 

In Appendix A, we provide the detection limits computed with 
the rms-based method and the LPA ones, when computed. The 
corresponding LPA-derived values for periods of 3, 10, 100, 500 
and 1000 days are reported in Table 2. We recall that the LPA 
method provides more accurate and better detection limits than 
the rms-based method. The rms-based limits are provided here 
only because for several stars, we do not have enough data to 
compute the LPA detection limits. In Appendix A, we also in- 
dicate the rms-based detection limits that would theoretically 
be achieved with a much better temporal sampling ("achievable 
rms-based detection limits"), as well as the rms-based detection 
limits that would be achieved if the stars were not active or pul- 
sating, i.e. if the noise was dominated by photon/instrumental 
noises. 

The detection limits depend on both the monitoring qual- 
ity (time-span and sampling) and on the stars (spectral-type, 
v sin/ and activity/pulsations). For most stars, the time-spans 
of the RV monitoring allow us to compute significant detec- 
tion limits for periods smaller than a few hundreds days. For 
a few targets, the time-span is close to or larger than 1000 
days. This is the case for HD39060, HD71555, HD105690, 
HD146624, HD172555, HD174429, HD188228, HD216956, 
and HD2 18396. 

Globally, the conservative rms-based achievable limits are 
all below 10 Mj up for periods less than 100 days, except for the 
A-type highly pulsating stars HD71555 and HD218396. As ex- 
pected, the LPA method provides detection limits better than the 
rms-based limits for all our stars. 

For early-type stars, the detection limits strongly depend on 
the available data. This is very clear if one compares HD39060 
(/JPictoris) and HD172555, which have similar spectral types 
and RV jitters. The detection limits can be substantially im- 
proved when averaging RV obtained on long series of data (1- 
2 hours). Planets of 2-3 Mj up can be found at periods up to 
100 days on HD39060 (see an extensive discussion in Lagrange 
et al., 2012), to be compared to masses of 4-10 M Jup for sim- 
ilar periods for HD172555. HD109536, with a temporal sam- 
pling covering 93 days only, has detection limits in the range 1-2 



Mj up for periods up to tens of days, even though its rms(RV) is 
as large as 216 m/s. We conclude that provided enough and well 
sampled data are available, it is possible to detect giant planets 
down to a few Jupiter masses, even with long periods, around 
these young early-type stars. 

Note that excellent detection limits are obtained for 
HD216956 (Fomalhaut) for periods up to more than 1000 days 
(Fig. 7). A ^ 30 Mj up brown dwarf at 5 AU was proposed by 
Chiang et al. (2009) to explain the observed Hipparcos acceler- 
ation. The present detection limits do not allow to rule out such 
a companion, because the time-span is not long enough, but the 
presence of such a companion could easily be tested provided 
a longer time-span and an appropriate observing strategy, as we 
did in the case of ySPictoris (Lagrange et al., 2012). 

Detection limits for late-type stars for periods below 100 
days are most of the time lower than 1 M Jup . Only for the few 
stars with the highest RV jitters, for example for HD 14 1593 
which has a RV jitter of 235 m/s, these detection limits become 
of the order of 5.5 M Jup for periods of 100 days. Stars with inter- 
mediate RV jitters, which are the most representative of our sam- 
ple, have much lower detection limits. For the well monitored 
HD105690, with a RV rms of 60 m/s, we obtain a detection limit 
of 1.2 Mj up for a period of 100 days, and of 2 Mj up for periods 
around 800 days. HD 183414, with RV rms of 60 m/s, has sub- 
Jupiter mass limits up to 100 days. Finally, HD181327, which 
has a lower RV jitter of 15 m/s, has a detection limit between 
less than 0.1 Mj up and 0.2 Mj up for periods up to 40 days, and 
0.3 M Jup for a period of 100 days. Note that the periods consid- 
ered here are limited only by the time spans of the data available. 
Additional data would allow to increase the range of periods to 
be considered. 

We conclude that when enough and well-sampled data are 
available, it is possible to find giant planets down to very low 
masses, down to 1 M Jup or below, orbiting these active young 
late-type stars, up to long periods. 

4.3.2. Improvements of the detection limits 

As mentioned above, adopting adequate observing and averag- 
ing strategies is particularly efficient to improve the detection 
limits for early-type, pulsating stars. This is also illustrated in 
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Fig. 7 in the case of HD2 16956, for which we recorded at each 
epoch long data sets of 1-2 hours. In this case, the RV jitter is 
reduced from 56 m/s to 25 m/s after averaging. 

For active stars, the detection limits obtained when taking 
into account the temporal structure of the stellar noise are al- 
ready very satisfactory. Nonetheless, it is worth trying to im- 
prove our detection capabilities. Several approaches have been 
proposed and gave more or less successful results on solar-type 
MS stars: use of simultaneous photometric variations to correct 
from spots/plages (Aigrain et al., 2012, Lanza et al., 2011, and 
references therein), use of the BVS (Boisse et al., 2011, and 
references therein) and use of spectroscopic activity indicators 
(Boisse et al., 2011). Using simultaneous photometric data may 
be quite powerful when the signal is not dominated by convec- 
tion (see an analysis in Meunier & Lagrange, 2013). Using the 
BVS or other spectroscopic diagnostics is very tempting as these 
parameters are obtained exactly simulatenously with the RV data 
and do not request additional observations. 

As seen previously, we do not observe any correlation be- 
tween the Ca emission and the RV variations, so we cannot use 
the Ca emissions to improve our detection limits. We then con- 
sider the use of the BVS. The v sin; of our stars are usu ally sig- 
nificantly higher than the spectral resolution, and this induces 
complex relationships between the BVS and the RV, as demon- 
strated by Desort et al. (2007) and Boisse et al. (2011). This is 
also seen in Appendix B in a number of examples. Departures 
from a linear correlation such as inclined "eightshapes" or even 
more complex structures can be observed. The amount of depar- 
ture from a simple linear dependance depends on the star v sin; , 
on the star orientation with respect to the line of sight sin ;' and 
on the spots parameters. 

We focused on the three active stars with the largest num- 
ber of data and for which a good (RV, BVS) correlation trend 
(hereafter slope) could be determined, namely (with the mea- 
sured RV jitter and v sin; in parenthesis): HD141943 (235 m/s, 
35 km/s), HD181327 (15 m/s, 17.9 km/s) and HD183414 (67 
m/s, 11.1 km/s). The slopes of the (RV, BVS) correlation are in- 
dicated in Table 2. After correction of the RV using the (BVS, 
RV) correlation, the RV jitters are significantly reduced, as well 
as the detection limits (see Fig. 8). For HD 14 1943, the RV jitter 
of the corrected RV is reduced from 238 m/s to 131 m/s, and 
the detection limits are reduced by a factor of 1.6 on average; 
they now fall in the range 1-5 Mj up for periods below 100 days. 
For HD183414, the RV jitter is reduced from 67 m/s to 25 m/s, 
and the detection limits are reduced by a factor larger than 2.9 
on average; they are at the Saturn mass level for periods below 
100 days. Finally, for HD181327, the RV jitter is reduced from 
15 m/s to 12 m/s, and the detection limits are reduced by a fac- 
tor of 1.5; they are in the 2-4 Neptune mass domain for periods 
below 100 days. 

An even better correction could be expected by removing 
the spot signal itself, using for example simultaneous photomet- 
ric data. A perfect correction would remove all activity-induced 
RV signal, and the achievable detections limits would be those 
achievable under photon and instrumental noises only (see the 
third panels of Fig. A.l). As an exercice and for illustration pur- 
poses only, we tried to fit the observed RV of the same stars 
with a keplerian signal (mimicking a spot) and then estimate the 
detection limits on the residuals. Significant improvements were 
obtained for HD141943 and HD181327, as seen in Fig. 8: the de- 
tection limits after correction now reach levels below 1 Mj up for 
the very active star HD141943 and of 1-2 Neptune mass planets 
for the weakly active star HD181327. However, this approach 
assumes that the spot signal can be fitted by a keplerian signal, 



which is not always the case. In particular, this assumption does 
not apply when the spot is not always visible or when the activ- 
ity pattern is complex, as illustrated by the simulations showed 
above. Hence, even though the detection limits are significantly 
improved, we would not advise to use this method for planet 
detection unless simultaneous photometric data are available to 
provide independant constraints on the spot(s) properties. 

4.3.3. Comparison with previous surveys 

Before comparing our detection limits to those resulting from 
other works, we recall that these quantities strongly depend on 
the targets themselves (jitter and temporal distribution of the 
stellar noise), on the observing strategy, and also on the way 
the limits are measured, as extensively discussed in Meunier 
et al. (2012). Our detection limits are much better than those 
observed by Bailey et al. (2012). However, the later focused on 
stars with spectral types M, i.e. later than those of our targets, 
and their detection limit estimation is probably conservative (like 
our rms-based method). Our detection limits, even before cor- 
recting from activity, are better than those obtained by Paulson & 
Yelda (2006) on the [3 Pic moving group; for these targets, their 
detection limits at short periods fall in the brown dwarf regime, 
where as our limits are well in the planet regime for similar peri- 
ods. We attribute this to the internal precision of their measure- 
ments, their adopted observing strategy, and their approach to 
estimate the detection limits. Finally, the fact that we find much 
lower detection limits for older stars than for younger ones is in 
agreement with Paulson et al. (2004) results. 

4.4. Search for circumstellar Ca in rapid rotators 

Finally, in no cases did we find evidence of circumstellar lines, 
except in the case of /JPictoris. In the case of the A-type 
ySPictoris analogous star HD172555, a narrow and faint compo- 
nent is present at the bottom of the Ca line. However, it is most 
probably of interstellar origine. 



5. Towards a full exploration of young stars giant 
planet population 

We have shown that RV monitoring with an appropriate observ- 
ing strategy and data analysis can allow to detect giant planets 
around young nearby stars despite their activity. The separations 
of detectable planets with RV techniques are however limited 
to a few AU. To get a complete view of planetary systems, it 
is necessary to complement RV data with other types of data: 
astrometric ones, e.g. Gaia, for typical separations of 2-5 AU 
(Sozzetti, 2011), and, above all, direct imaging data which can 
provide detection limits at much larger separations. Direct imag- 
ing is particularly interesting for young systems as the contrasts 
between young planets and their parent stars are more favor- 
able than for older systems. The combination of these techniques 
therefore offers a unique opportunity to explore the giant planet 
population on a same class of stars. We illustrate in this section 
the potential of the coupling of RV and AO deep imaging data, 
using the example of Harps and NaCo, and Sphere data. 

We consider three stars of our sample for which we have 
both Harps RV and NaCo imaging data, and which could 
be considered as prototypes for such complementary studies, 
namely HD183414, HD181327, and HD188228. The G-type 
star HD183414 is an example of medium activity-level star, 
with a jitter of 67m/s; it is located 35 pc away. The F-type star 
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Fig. 7. Impact of data averaging on pulsating stars: rms-based detection limits obtained on HD 216956 when using all measurements directly (left 
panel), and when averaging the data over 2 days (right panel). Same symbols as for the middle panels of Fig.A.l. 



HD181327, located at 51 pc, has a low activity level (jitter=15 
m/s). The pulsating A-type star HD 188228, located at 32 pc, has 
a jitter of 127 m/s (amplitude of 2 km/s). We plot in Fig. 9 the 
RV and imaging detection limits expressed in Jupiter masses for 
these three stars. The deep imaging detection limits are either 
those measured on NaCo data (Rameau et al., 2013), and the 
ones expected with Sphere (Mesa et al., 201 1). Note that the de- 
tection limits of HD 183414 are those obtained once correction 
for activity using the RV-BVS relation has been made. Different 
star inclinations are considered for the RV data. The limits are 
strongly impacted for inclinations below 60 degrees. For larger 
inclinations, the RV derived detection limits are very similar. 
We see that the detection limits are in the planetary regime with 
both approaches. In these examples, the RV limits are much bet- 
ter than the NaCo ones for the G-type stars (this effect is even 
stronger when considering moderately active stars with a RV jit- 
ter below 50 m/s), while the AO limits are better than the RV 
ones for the A-type star. 



The complementarity of HARPS and NaCo in terms of sepa- 
ration ranges is obvious. However, gaps are still present. This is 
particularly true for HD 183414, for which the available RV data 
are spread over 189 days only. On the contrary, for HD 188288, 
which is located at a similar distance, there is an excellent over- 
lap because we dispose of RV data spread over about 2000 days. 
Finally, the gap observed for HD 18 1327 is smaller than that of 
HD183414 because the available RV data span is much larger 
(560 days), even though the star is located further away. We note 
the important impact of forthcoming imagers like SPHERE on 
the VLT: with inner working angles smaller than 0. 15"-0.2", they 
will allow the detection of planets much closer to the stars than 
current imagers such as NaCo; this will fill the gap between RV 
data and imaging ones for stars closer than typically 30 pc; this 
will also significantly contribute to reduce the gap for stars lo- 
cated further away. Finally, we see that Gaia can help filling the 



gap, if any, between RV and AO direct imaging data. This will 
be particularly important for distant stars. 



6. Concluding remarks 

We have shown that RV data allow the detection of giant planets 
around stars in close-by, young associations. Given the available 
data, we could explore different separations ranges, up to more 
than 2 AU for some stars of our sample. For low to moderately 
active stars, detection limits of a few Neptune masses are ob- 
tained. The best detection limits are obtained for low to moder- 
ate v sin/ in the range 10-30 km/s, which are generally observed 
for stars older than 30 Myr, although there are some exceptions. 
Early-type stars are generally pulsating, and averaging over the 
high frequency variations reduces considerably the achievable 
detection limits. Values well within the giant planet domain are 
obtained. 

Finally, we have illustrated the tremendous interest of cou- 
pling RV and deep AO data using Harps and NaCo data, and, 
even more, future Sphere data. RV data covering more than 2000 
days allow planet detection and characterization for separations 
up to 3-5 AU, while forthcoming imagers will provide excellent 
detection limits for separations greater than 0.15-0.2", which 
correspond to 4-5 AU at 30 pc. For more distant stars, Gaia will 
fill the remaining gap. We will have in the next decade, and for 
the first time, a unique opportunity to make an exhaustive giant 
planet search around these stars, that will in turn lead to unprece- 
dented studies of giant planet formation and evolution. 
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Fig. 9. Upper panel: AO detection limits (observed NaCo limits in plain lines and expected SPHERE limits in dotted lines) and Harps RV detection 
limits for the active G-type star HD183414 (192 Myr, 35 pc and time-span of 189 days). In both cases, different star inclinations: 10 (light blue), 30 
(red), 60 (blue) and 80 degrees (green) are considered. The Gaia detection limits are also indicated. Middle panel: idem for the less active F6 star 
HD181327 (12 Myr, 51 pc and time-span of 588 days). Lower panel: same for the pulsating A-type star HD 188228 (40 Myr, 32 pc and time-span 
of 1984 days). Note that as the star is brighter than V=6, no Gaia data will be available for this star. 
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Fig. 8. Upper panel: LPA detection limits after correction of the RV 
using the (BVS,RV) slope (blue dots), and after fitting and removing 
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35 km/s). The uncorrected LPA detection limits are indicated by green 
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Appendix A: Observed RV, BVS and measured 
detection limits 

Appendix B: Simulations results 
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[i-siu] A>I 



[,-SI0] A« 



[,-Sin] Atf 



[i-sui] A>I 



A.-M. Lagrange et al.: Planets around stars in young nearby associations, Online Material p 7 




[%V] (Ows-jv 




io o <n 



[,s ui] ireds jopsstg 



I i i i i i i i i i | i i i i i i i i i | 



I: 

C-; _ 

I - 



o 



- CM 



O 



A.-M. Lagrange et al.: Planets around stars in young nearby associations, Online Material p 8 





RV 


[m s 1 ] = f(JD-2453000 [days]) 


4(1 




i i i | i i i i i^i i i | i j 


20 




.* • " 





-.... 





20 


- 


/2 


-40 


-I i 


, i , ,;•-,-; 



Bisector norm flux = f(RV [m s 1 ]) 





3isector span [m s" 1 ] = f(RV [m s" 1 ]) 



1000 1005 1010 1015 



480 500 520 540 




Bisector span [m s 1 ] = f(JD-2453000) 




1000 1005 1010 1015 





RV [m s- ! 


= f(JD-2453000 [days]) 








00 




-_ 













- 








00 




„ 1 






« . 




"i , i 


1 , 1 , 1 ** , I - 



Bisector norm flux = f(RV [m s 1 ]) 



Bisector span [m s 1 ] = f(RV [m s 1 ]) 






Bisector 


span [m s 1 ] = f(JD-2453000) 


100 



100 




i ' | ' i ' L, ' i- 



RV [m s'] = f(JD-2453000 [days]) 



Bisector norm flux = f(RV [m s" 1 ]) 



Bisector span [m s" 1 ] = f(RV [m s" 1 ]) 





1.0000 1.0005 1.0010 1.0015 
10+ 3 





RV [m s 


] = f(JD-2453000[days]) 






i ' ' '«w I_ 








500 




- 














....... • . - 






. _ 


500 




• * : 






. . 




"1 , i , 


1 «f T , 1" 



1000 2000 

Bisector norm flux = f(RV [m s -1 ]) 




-1000 1000 

Bisector span [m s~'] = f(RV [m s 1 ]) 





3isector 


span 


ms" 1 


= f(JD-2453000) 


500 



















.*** 








-,«*- 


500 








i 1 r% i i 1 i T 



1.0000 1.0005 1.0010 1.0015 
10 +3 



Bisector span [m s'] = f(JD-2453000) 



0.98 


11111111 ''StV^CWsk 1 ' - 


0.96 




'M '■ 


0.94 




'11 1 


A\k ' 


0.92 




iff 


\/NX '■ 


0.90 


r^CX-r^K\ I 1 1 


ElfoCKi 



1000 




i i j.i i i i i i i i i 


-0.8 : 


500 



: 


*. 
n 

\ 

--X-. 


-_ 


-500 


- 




- 


1000 


r, , 


"CI 


_i_ji 



1.0000 



1000 2000 3000 






RV [m s 1 ] 


= f(JD-2453000 


[days]) 










1000 




* 


* 




- 


. 


* 




■ 


• 


- 













"**"*""* 


*********** 


* .~~ 




■ 




• - ~- 










-1000 


- 




• . 




i , i 


i.i.i 


i **. i 





Bisector norm flux = f(RV[m s" 1 ]) 




0.98 




- 


0.96 




-_ 


0.94 




) ~_ 


0.92 


f i id i/ 1 iSll fi i\ *U i\ i ><1 i\i \ 


~>T 





Bisector span [m s" 1 ] = f(RV [m s" 1 ]) 


1000 


-' 


.1 .| 1 1 1 1 | 1 1 1 1 | 1 1 1 L 






-0.5 - 








500 




^ . - 






* ^ . 













- .*" 












^ 


-500 




" 






. . - 


-1000 




- 









1.0005 1.0010 

10+ 3 

Bisector span [m s'] = f(JD-2453000) 




1.0000 1.0002 1.0004 1.0006 

10« 



1000 2000 3000 4000 



1.0000 1.0002 1.0004 1.0006 

10 +3 



Fig. B.l. First column: RV versus time of a simulated equatorial spot covering 1.5% of the visible surface and a temperature contrast of 1200 K, 
for a seen edge-on sun-like star and for five different values of v sin; from the upper panels to the lower panels: 3 km/s, 10 km/s, 30 km/s, 50 km/s, 
80 km/s. Second column: same for the bissectors. Third column: BVS versus RV. Fourth column: BVS versus time. 
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Fig. B.2. Same as Fig. 10 for an inclination of 30 degrees (three first lines) and 5 degrees (three last lines), a colatitude (from the pole) of 30 
degrees (except 5th line with a colatitude of 60 degrees) and v sin/ (from upper panels to lower panels) of 3, 10, 40, 10, 10 and 10 km/s. 



